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ABSTRACT 

Mellen  and  Moffett  (-1,  3)  have  derived  models  for  parametric  sonar 
radiator  design  following  Westervelt's  (2)  theory.  Their  models  analyze 
both  absorption  and  saturation  limits  in  both  the  near  field  and  far  field. 
This  memorandum  is  concerned  with  the  far  field  zone  expansion  of  the  far 
field  beam  pattern  and  pressure  into  a  series  of  eigenfunctions  which  can 
be  integrated  term  by  term.  After  summing  over  the  terms,  we-can  show  that 
the  far  field  beam  pattern  is  a  product  of  the  two  primary  beam  patterns. 
The  conversion  efficiency  reduces  to  an  experimental  integral  which  appears 
in  Reference  (b). 


ADMINISTRATIVE  INFORMATION 

This  memorandum  was  prepared  under  Project  No.  A65405,  Task  No. 
ZR00Q101,  "Parametric  Sonar  Research,"  Principal  Investigator  Dr.  R.  H. 
Mellon,  rprie  PA41 ,  The  sponsoring  activity  is  CNM/MAT  C33,  Dr.  J.  H. 
Hath. 


The  author  of  this  memorandum  is  located  at  the  Newport  Laboratory, 
Naval  Underwater  Systems  Center,  Newport,  Rhode  Uiand  02840. 
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INTRODUCTION 


The  purpose  of  this  memo  is  to  present  in  seme  detail  a  derivation 
of  the  far  field  zone  beam  pattern  of  a  parametric  sonar.  Moffett*  has 
derived  from  Westervelt's  theory^  the  two  fold  expression 
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for  the  pressure  of  the  secondary  v/ave  where  R  is  the  primary  wave 
collimation  distance.  The  volume  integral  is  to  be  performed  over  the 
region  occupied  by  the  primary  wave.  In  this  expression  £  is  the  distance 
(r^rO  between  the  field  point  and  the  integration  point.  The  other  terms 
are  defined  in  the  section  "Glossary  of  Terms." 

The  first  integral  in  Equation  (l)  represents  the  contribution  to  the 
far  field  from  near  field  interactions.  It  is  evaluated  in  references  (1) 
and  (3).  The  essence  of  these  papers  is  that  the  secondary  beam  pattern  is 
a  product  of  the  two  primary  beam  patterns.  The  parametric  efficiency!  rP  I 
is  derived  in  both  references  and  the  directivity  index  DI  is  derived  R0P0 * 
in  reference  (3).  Curves  of  the  directivity  index  end  parametric  efficiency 
are  presented  in  reference  (3). 

The  second  integral  represents  the  contribution  to  the  far  field  from 
the  spherically  spreading  region.  It  is  the  intent  of  this  memorandum  to  . 
evaluate  this  integral. 
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FAR  FIELD  INTERACTION  ZONE 

We  rewrite  the  second  integral  in  Equation  (l)  in  the  form 


0Po2R2k2eiwt 
P(r)  -  _ 
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where  4  has  been  replaced  by  jr-r^ji^e  distance  between  the  field  point  r 
and  the  integration  point  t'.  S2(r')  is  the  source  function.  The  volume 

element  dv'  is  r'2  sin  dr'  dO*  d$'and  r'  extends  from  R2  to®.  The  distance 
R2  is  defined  to  be  k0R0  and  is  called  the  spherical  wave  distance. 

“IT" 

The  justification  for  extending  the  range  of  integration  from  R  to  R2 
was  that  most  of  the  contribution  in  the  first  integral  in  Equation  (1)  came 
from  on  the  axis  and  the  integrand  was  non  oscillatory  up  to  r'=  R2»  Thus  the 
near  field  contribution  was  extended  out  to  the  spiericai  wave  distance  ana 
the  far  field  or  spherically  spreading  form  was  assumed  to  begin  at  P.2. 

1 _ 

The  quantity  exp  {-ik|r-ri  )  is  the  Green's  function  for  the  Helmholtz 
operator.  It  has  the  expansion  1 

e~ik (r-r  )  _ji.  9  ro-.. n  n  t _ «•»  •  /1  .  f ..  .. 
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where  Pn(co$  0)  »  PB(cos0  cos  0"  +sin  0  sin  0'  cos  if') 


P„(cos  0)  P_ (cos  0')  ♦  21  (n-m)J 

n  maiTKTjiit: 
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The  terms  Pnm(cos 9 )  are  the  asssociated  Legendre  polynomials  of  order  n. 
If  we  integrate  over  the  azimuthal  angle  <j>  ,  we  have  from  (4) 


J  PQ(cos  0)  •  2irPn(cos  0)  Pn(cos0"). 


In  the  far  field  we  use  the  second  of  the  expansions  in  Equation  (3). 
Combining  the  results  of  Equation  (5)  with  the  expansion  in  Equation  (3) 
we  rewrite  Equation  (2)  as 


r  P  «  BP02k2R02eiwt 


j  j* (-ik)  2  (2n+l)  P  (cos  0)  P  (cos  0")  x 

J  J  n-0  n  nv 


Mkr')  hn(2)(kr)  e"ikr*  JP(0')|*sin  e'  dr' 


de". 


Thp  source  function  S2(r")  has  been  resolved  into  a  product  T^(r ')  )|  2 

where  T^(r')  is  a  linear  taper  function  and  |  F(  Gf) p  is  the  product  of  che  two 
primary  beam  patterns.  In  the  far  fj.eld  Hankel  function  hn'2'(kr)  can  be 


written  as 


h  (2^fkr}  =  e“i  r<r”£(n+lH]  =  in+1  e”ikr 

n  —  17~ 


We  then  have 
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This  expression  con  be  written  symbolically  as 
-ikr 

n 
kr 


Ke 


Z  (2n*l)  in  anbn  Pn  (cos  0) . 
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where 
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(9) 

(10) 


The  inlegiand  has  been  resolved  into  a  series  of  eigenfunctions 
which  can  be  integrated  term  by  term.  We  proceed  first  with  the 
integration  over  the  angle  0  . 


Each  of  the  primary  beams  is  assumed  to  have  the  beam  pattern 


2J , (k  a  sin0) 
F(0")  «  1  0 


(H) 


kQ  a  sin  0" 

where  a  is  the  radius  of  the  transducer  and  fcQ  is  the  average 
primary  wave  number.  The  aperture  a  is  related  to  the  primary 
wave  collimation  distance  RQ  by  , 
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for  a  circular  pistons  the  coefficients  b  Becomes  in  Equation 

00) 


b 

n 


*  IVk0a  sin  6'j 

L  k.a  sin  0'  _ 


sin  6  dQ'. 


(13) 


The  square  of  the  Bessel  function  has  the  representation 


*  -  (-l)n(2m*2)‘  (.Vx)Im+2 

J,  (x)  -  l  _ 

m-0 

n «  (m+2) l  [ (m*l) !] 2 


Thus  b  becomes 
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Let  x  =  cos  0  .  Then 
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In  this  case  HQ(2koa)  and  H|(2k0a)  are -Struve  functions  of 
orders  0  and  I  and  argument  2k  a.  Using  the  Lagronge  Duplica¬ 
tion  formula,  we  can  representt^n  as  follows: 


(-l)n  2i_2n  (2n)«  r(3/2) 


(n!)  r(l-n)  r(n+3/2) 


3  1,  1;  3,  2,  n+3/2,  1-n;  -4k  2 


where  3^4(0}/  <*2>  °3/  P]/  P 2 •  P3>  P4J  z)  *s  a  generalized 
hypergeometric  function. 


Returning  to  Equation  (9)  we  calculate 


00 

-an  -  f  T2<r')  jn  (hr')  e"lkr' dr', 

•'r 

2 

where  n  is  even  since  b2  +|  =  0.  Following  the  work  of  Mellen 
and  Moffett  ,  v/e  choose 


T  (r)  = 


1+(X  sinh  ~1  r'/Ro)2 


••  rv  r  * 

.The  term  e“  represents  the  abscrbtion  of  lie  primary  wave 
and  the  term  tl  +  (X/2sinh  rVR0)2]  -1/2  js  due  to  the  propa¬ 

gation  of  the  saw  tooth  .wave ,  (see  Equation  (;9)  of  reference 

(I)). 
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If  we  let  x'  =  R^P.  we  have 
o 


r  e~2a  V 

1/ai+x  sinh“lp 


e“iokoRop  j2n(akoR  p)  dp, 


where  A  represents  the  downshaft  ratio  fQ  =  k  .  The  first 


two  spherical  Bessel  functions  are 


T*  TT 


Jo(z)  = 


sin  z 


z  3  ok  R  p 
o  o 


j2(0 


(3  -  _1)  sin  z  -  3_ 
z7  z  I2 


cos  z 


~  __  sin  z  +  .3  (stn  z  r  cos 
z  ~z3 


It  is  expected  that  because  the  variable  of  integration  is  iarge 
to  start  with  that  the  term  in  parenthesis  in  |2  will  contribute 
very  little  to  the  integral  since  each  of  the  terms  is  oscillatory 
and  cancels  each  other  out.  It  is  expected  that  |2n(z)  can  be 
set  approximately  to  (-l)n  jc(z). 

When  the  saturation  parameter  x  In  Equation  (22)  is  non 
zero,  the  integral  cannot  be  evaluated  in  closed  form.  A 
series  of  machine  computations  was  made  by  Mr.  Marvin  Goldstein 
of  the  New  London  computer  laboratory  for  tho  following  sets 
of  parameters  as  shown  in  Table  I. 


TM  No. 
PA4-3  13-71 


TABLE  I 


Parameters 

Used  in 

E  va  1  uati  n 

X 

1 

1 

a 

2aRc 

0 

5 

1 

CM 

O 

• 

10 

i0 

.04 

20 

100 

.06 

50 

1000 

• 

o 

CO 

.  10 

s  case  1  /'o  = 

10  and 

1  =  100 

M 

o  o 


100 


shown  in  Table  II. 


2a  R, 


=  100,  the  following  results  arc 


TABLE  II 
Values  of  02n 


2 

2 

2 

2 

X 

ao 

a2 

a4 

a6 

0 

. 8271x10"2 

.  8271x10-2 

.  827  ixl0-2 

. 8256xl0"2 

02 

. 8240x20-2 

.  824  1 x  1 0 ”2 

.  824  IxlO-2 

. 8226xl0"2 

04 

.8150x10“2 

.  8 1 5 lx 10“2 

.8151x10-2 

. S136x  10"2 

06 

.  8004x 1 0“ 2 

.  8004x  10-2 

.  800''x  1 0 “ 2 

. 7990xl0-2 
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The  values  of  |  aQ  2#  |o2|^#  and  ja^  2  agree  to  1  part  in 
104  apd  the  value  of  ja^j  2  agrees  witjj  the  others  to  15  parts 
in  10  .  However,  the  series  for  |a  j  was  converging  so 
slowly  that  the  agreem''^'  could  be°close-  As  a  very  good 
approximation  (and  r*«-  which  was  expects  •'  all  values  of  | a 2 r  | ^ 
can  be  expected  to  be  the  same  as  |aj  2.  We  now  return  to 
Equation  (8),  The  fact  that  aj,  03,  ...  03^+]  was  not  calculated 
is  not  important  since  the  b*  ,,  is  identically  0.  We  can  v/rite 
Equation  (8)  as  * 


P  -  E  (2n+l)  in  (-l)n  b  P  (cos  0) 

kr  0  n  n  n 


since  a2n  =  (-l)n  a0.  This  becomes 
K  -ikr 

P  *  — —  2a  E  /  lF(0')|2  (2n+l)  P  (cos  e)  P  (cos')  sin  0'd9',  (24) 

-  Kr  0  n  Ja  0  n  n 


The  terms  (2n  + 1 )  Pn(cos  d)  Pn(cos  0-")  are  normalized  signfunctions 
and  by  the  closure  property 

nf0  un<x)  »*<x')  =  fi(x-x')  (2 


Thus 


V  -ikr 

P  *  M -  2a 

kr 


„/  lF<e')l: 


6  (cos  8'-  cos  0')  sin  8  '  d0  ' 


2a0|P(9)| 


’  )2 


TM  No. 
PA4-3 1 3-7 1 


Strictly  speaking,  the  beam  pattern  in  the  far  field  is  a 
very  weak  convolution  of  the  two  primary  beam  patterns  but 
to  a  very  good  approximation  as  shown  in  Equation  (27),  the 
beam  patterns  >  n  the  far  field  is  a  product  of  the  two  primary 
beam  patterns. 

If  one  considers  a  value  of  the  saturation  parameter  x  to  be 
zero.  Equation  (22)  can  be  integrated  in  closed  form.  Let  us 
write  the  pressure  distribution  as 


p  _  iwt  -ikr 

r  -  u  e  e 


pc 

io^<rkoR0P)  dP' 


kr 


|F(0)|2  kRc  j[7e"2aRoP‘'  Mo* 


arrd  consider 

»  CC 


f  #-2.R.p-l,R„p  jo(rk  Rop)  dp 

Micr 

=  f"  e'2aRo P  dp  -  C  e'2oR°p  "  2'crkr»Rr>P 

Jl/n-  - ~ -  •'1/tr 


dP» 


2i<rk  R  p 


2jaRfi  P 


2 1  cr  k  R 
o  o 


^  E'(^!l+  2ik“R| 


For  imaginary  k  R  =  2  00  ,  z  eZE|(z)  =■  I. 


o  o 


Thus  E1/2a^H2ikoRoWex{^i°"2,koRo  j 
'  a  '  2aR_  «...  ■“ 


o  -2i k  R 
—  oo 
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Hence 
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P  =  PP°2|i02|<2  eiwt  e-,kr  !f(0)  2  kR  1 

o  k  r  I  ° 


E,(2aRo) 


2 !  o-  k  R 
o  o 


Taking  the  value  of  the  parametric  efficiency  r P  ,  we  have 


RoP0' 


1  R  p 

or  o 


E,(2aRo)  j  F(6)| 


which  is  the  same  value  as  Equation  35  of  Reference  I.  Equation 
(35)  of  Reference  I  is  the  on  axis  value  of  the  far  field  parametric 
efficiency. 

CONCLUSIONS 


This  memorandum  has  shown  that  the  secondary  beam  pattern 
of  the  far  field  interaction  zone  is  a  very  week  convolution  of 
the  two  primary  beam  patterns  which  to  a  very  good  approx i ma t i on 
is  a  product  of  the  two  primary  beam  patterns.  The  pressure  ampli¬ 
tude  of  the  secondary  wave  was  calculated  using  all  the  off  axis 
contributions  and  agrees  with  that  calculation  made  on  the  basis 
of  the  on  axis  contributions  in  References  I  and  3.  This  gives 
great  credence  to  the  approximations  made  in  these  two  references 
so  that  all  the  quantities  and  curves  appearing  in  these  two 
references  can  be  considered  to  have  much  validity. 
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l  I 

r-r'j  =  Distance  between  field  point  and 
integration  poi nt 


Densi  t  y  of  water 


Reciprocal 


of  downshift  ratio  = 


Saturation  number  = 


PPokoR 


o 


Angular  frequency  of  secondary  =  2nf 
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